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Chapitre 4 

Interpolation et approximation 
polynômiale 

4. 1 Introduction 

Tout au long des deux chapitres précédents, nous nous sommes intéressé principalement au problème 
de l'approximation d'un nombre (zéro d'un polynôme ou le zéro d'une fonction non linéaire). Nous 
abordons dans ce chapitre un nouveau type de problème faisant intervenir la notion d'approximation 
d'une fonction. 

La méthode usuelle d'approximation des fonctions est la méthode d'approximation polynômiale. 
Parmis les divers types d'approximations polynômiales utilisées, la plus flexible et la plus facile à 
construire est l'approximation par interpolation polynômiale. 

4.1.1 Objectif de l'interpolation 

Etant donnée une fonction f(x) = e x (ou autre expression de / beaucoup plus compliquée) que 
nous désirons approcher par un polynôme de degré n = 2 (ou n > 2)) sur l'intervalle [—1, 1]. Nous 
allons voir que nous pourrons construire trois polynômes à ce sujet: 

• Le polynôme de Taylor pt{ x ) — \ x2 + x + 1- L'erreur maximale de cette approximation est 
0.21828. 

• Le polynôme de Lagrange passant par (— l,e ), (0,1) et (l,e) est le polynôme de dgré 2 
donné par: pl( x ) — 0.54308a: 2 + 1.1752a: + 1. L'erreur maximale de cette approximation est 

0.174378. 

• Le polynôme de Newton passant par ( — 1, e ), (0, 1) et (1, e) est: 

p N {x) = 0.367879 + 0.632121(.t + 1) + 0.543081(a; + l)(x - 0). 

L'erreur maximale de cette approximation est 0.174378. 

Dans tous les calculs portant sur / (tels la dérivée de / ou l'intégrale de /) on pourra remplacer 
/ par Pt ou Pl (ou par un autre polynôme). En fait il y a plusieurs méthodes pour construire un 
tel polynôme. Dans un premier temps nous allons nous focaliser sur l' approximation de TAYLOR, 
Lagrange, Newton et celle de Chebyshev. 
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4.1.2 Un exemple constructif 

Un résultat classique en analyse est le développement en série entière des fonctions usuelles telles 
que sin, cos, ... 



c (2n+l) 

(2n + l) 

E (2n) 

<2n) 



-^ = E(-i) n kôï = J - - - 



X X 

2!~ + 4T <;! 

2 „3 ^,4 



ïxix = > ( — l) 71 - — = a; -H ; r + ••• pour tout x 6 IR 

f^ v ; (2n + l)} 3! 5! 7! A 

pour tout x € IR. 
u: 

e^V— - = 1 + 2; + — + — + — + ... pour tout x G IR 

*— ' n! 2! 3! 4! 

71 = 

234 

ln(l + #)=:£ (- — h . . . pour tout x g] — 1, 1] 

2 3 4 

Notre but est de savoir, jusqu 'a quel ordre ri (fini) nous devons écrire le développement pour avoir 
une bonne approximation de la somme infinie. 

Pour illustration, soit à calculer e = exp(l) à partir de l'écriture de (e x ) x= i en série entière. On 
a donc 

111 1 

ÏÏ + - 



: *-L n\ l+ 1! ' 2! ' 3! 

Tl— 



Si nous tronçons la série (infinie) précédente à un certain ordre n > 1, nous obtenons une approxi- 
mation de la valeur de e par: 

111 1 

5 " = 1 + Ï! + 2^ + 3! + --- + ^ 



n 


s„ 





1 


1 


2.0 


2 


2.5 


3 


2.666666666666 


4 


2.708333333333 


5 


2.716666666666 


6 


2.718055555555 


7 


2.718253968254 


8 


2.718278769841 


9 


2.718281525573 


10 


2.718281801146 


11 


2.718281826199 


12 


2.718281828286 


13 


2.718281828447 


14 


2.718281828458 


15 


2.718281828459 



La valeur numérique de e, à une précision de HT 29 , est donnée par e w 2.71828182845904523536028747135. 
On remarque donc, à partir du tableau ci dessus, que plus l'ordre n du développement est élevé 
plus l'approximation de e par la série de Taylor est meilleure. Sig donne une approximation de e 
avec une précision de 10~ 13 . 
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4.2 Approximation polynômiale de Taylor 

Théorème: 

Soient a et h deux nombres réels tels que a < b et / : [a, b] — > IR une fonction n fois dérivable sur 
[a, 6]. On suppose que / ÉE C' n+1 '[a, 6] et xo G [a, b] (un point fixe). Si x G [a,fr], alors il existe un 
élément c entre xq et :c tel que 

f(x)=P n (x) + E n (x) 



-PrcW = /(ïo) + 1! [x - xq) + ■■■ H -j (s - x ) 

f(»+l)M 
^W = J (n + 1 )/ (^-^)" +1 (4.1) 

P n est le polynôme qui peut être utilisé comme approximation de /, E n {x) est l'erreur commise 
lors de cette approximation. 

Preuve: 

On pose: 

P n {x) = f(x ) + ^-(x -x ) + ... + ^-^-{x - x ) n 
1! n! 

P n est un polynôme de degré n et est indéfiniment dérivable. Considérons la fonction g définie par: 
g(t) = f(t) — P n (t) + X(t — Xo)"' , A est une constante telle que g(x) = 0. On a donc: 

g(t) = fit) - p n (t) + f" (a;) " ( Sj (t - *o) ( " +1) . 
{x — XQ)y n+i -> 

On vérifie facilement que: 

p»(*o) = f(x ) , K(x ) = f(x ) , p;;(x ) = f(x ) , ... , j*>)(xo) = f {n \x ) 

Pn est identiquement nulle, donc g est n-fois dérivable sur [a, b] et g( n > est dérivable sur ]a, 6[. 

De plus on a: 

g(xo) = g'(x ) = g"(x ) = ... = W(xo) = 

Alors on a g(xo) = g(x) = d'après le théorème de Rolle il existe c\ entre xq et x tel que g'(ci) = 0. 
On a alors g'{x^) — 9'( c i) — 0; donc il existe un élément C2 entre a:p et c\ tel que g"{c2) = 0. Par 
récurrence on construit donc une suite d'éléments {01)^^2 ... n de ]a.b[ tels que: 

a < c n < c n -i < . . . < c 2 < ci < b , g^(Ci) = , si i = 1, 2, . . . ,n 

On a alors g^ n '(xo) = g^ n '(c n ) = et d'après le théorème de Rolle il existe c entre xq et c n tel que 
ff (n+l)( c ) =0 . Or 

j(»+i)(c) = = /C»+l)(c) + f^'f^ + 1)! 

[X — XQ)V l+i > 

donc: f(x) — P n (x) = / +i}i ( x ~ x o) +1 • Ce qui démontre le théorème. 
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50 3. Interpolation et approximation polynômiale 

4.2.1 Exemple d'estimation de l'erreur 

Nous allons montré dans cet exemple pourquoi nous avons eu besoin de prendre n = 15 pour obtenir 
une approximation de e par S\§ avec une précision de 10~ 13 (voir tableau précédent). 
f(x) = e x , f( n '(x) = e x , le développement de Taylor au point xq = à l'ordre 15 de / est donné 
par f(x) = Pis(x) + E 15 (x). 

9 3 1 r i 

_ , . x z x ô x lù 

p wix) =l +x+ - + - + ... + - 

-Pis(l) — 5i5- L'erreur est donnée par: 

E li {x) = } y ,ce]0,l[ 

A partir de /^ 16 )(c) = e c et de 1 < e c < e < 3, on obtient une majoration de l'erreur (pour x=l): 
E 1 ^l) = i^M = ^< JL = 1.43384310- 

4.2.2 A propos de l'erreur 

Il est évident que: 

i la précision du polynôme de Taylor est d'autant meilleure que n est grand 
ii la précision de chaque polynôme de Taylor décroît lorsque x est loin de xq. 
Un compromis entre (i) et (ii) et de choisir: 

• n assez grand 

• restreindre le domaine de validiter de l'approximation de / par P n a un voisinage de xq 

Si on choisit x dans un intervalle de centre Xo et de rayon R (c.à.d [œ-ïCol < R), on obtient la 

majoration de l'erreur suivante: 

MR n+1 
\KrT G ur\ = \E n (x)\<j^— (4.2) 



AI <max{\f (n+1) (x)\ : x - R < x < x + R} 
On remarque alors que: 

• Si n est fixé, et \f^ n+1 \t)\ est bornée sur [xq — R, xq + R] : le majorant de l'erreur (4.2) est 
donc proportionnelle à R n+1 /{n + 1)! et décroit lorsque R — > 0. 

• Si R est fixé, le majorant de l'erreur (4.2) est inversement proportionnel à (n + 1)! et tend 
vers zéro lorsque n devient suffisament grand. Ceci est schématisé dans la figure suivante: 

Sur cette figure on remarque que l'on approche bien e x par P^ix) sur [—2,2], 
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5 

P21.X) 



-P4u) 
P3(x) 



P2(x) 



Figure 4.1: Approximation polynômiale de Taylor. 

4.3 Interpolation polynômiale de Lagrange 

Connaissant (n + 1) points du plan IR X 1R, Aq(xq, /o), Ai(xi, /i), ^2(^2, /2), • ■ ■ , A n (x n , f n ) tels 
que les Xi sont tous distincts. On cherche un polynôme P, de degré au plus n, passant par ces 
(n + 1) points c'est à dire: 

P{xo) = /o, P(xi) = /1, ■ ■ ■ , -P(x„) = /„ 

• Les fi sont souvent les valeurs d'une fonction / au points x%*. fi = f(xi) on dit alors que le 
polynôme P coincide avec la fonction / au point (#^=0,1,...™- Le polynôme P approche la 
fonction /, dans tous les calculs portant sur / on pourra la remplacer par P. 

• l'approximation de / par P est d'autant meilleure qu'au points (xi) nous aurons plus d'information 
sur / ainsi que sur /'. 



4.3.1 Existance et unicité du polynôme de Lagrange 

Théorème II existe un polynôme et un seul de degré au plus n tel que: 

P( Xi ) = fi , i = 0,l,...,n 

Preuve: En fait, il y a deux façons pour démontrer ce résultat: une directe et l'autre moyennat 
les polynômes de Lagrange. 

méthode directe: Soit P un polynôme de degré n ( au plus): 

P(x) = a n x n + a n _ix n ~ + . . . + a\x + ûq 

Par hypothèse on a: P{xi) = fi i = 0, 1, . . . , n, ce qui représente un système de n + 1 équations 
linéaires et n + 1 variables (les coefficients ao ai, . . . , a n de P), c'est un système de Kramer. La 
solubilité de ce système d'équation s'établit par l'étude de son déterminant: 

Xq Xq Xq 1 

xl x^' 1 ... Xl 1 



- ■ Xfl 1 



il s'agit du déterminant de Vandermonde qui vaut: 



^i=0,j>ii X i x j) 
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52 3. Interpolation et approximation polynômiale 

et qui est effectivement non nul du fait que les x% sont distincts deux à deux, d'où l'existance et 
l'unicité des a%. 

Formulation de Lagrange: 
Unicité: 

Si Pi et P2 sont deux polynômes de degré au plus n tels que: 

Fi(si) = P 2 (xi) = fi , i = 0,l,...,n 

Le polynôme D = P\ — P2 vérifie D(x%) = Pi(x{) — Pzix-i) = pour tout % = 0, 1, 2, .. . ,n. On con- 
clut que D de degré au plus égale à n s'annule n + 1 fois, D est donc le polynôme nulle donc Pi = P-2- 

Existance: 

Nous allons construire le polynôme de Lagrange qui satisfait aux conditions: 

P{x i ) = f i , i = 0,l,...,n 

questions préliminaires: 

• Q\. Trouver le polynôme L k de degré au plus n tel que 

L k (xi) = V i £ k (i = 0,1,.. . ,77.) L k (x k ) = 1 
L k s'annule n fois donc peut s'écrire comme: 

L k (X) = C(X - x )(X - ai) ... (X - Xfc-i)pC - sBt+i) . . . (X - z n ) 
ou C est une constante que l'on peut déterminer à partir de la contrainte L k (x k ) = 1. On trouve: 

C" 1 = (x k - x )(x k -xi)... (x k - x k _ 1 )(x k - Xk+i) ...{x k - x n ) 
finalement L k s'écrit comme: 

(X - x )(X - Xl )...(X- .t,_i)(X - x k+1 ) ...(X- xn) 



L k (X) : 



(x k - x )(x k -Xi)... (x k - Xk-l)(%k - x k+ i)... (x k - x n ) 



• Q2 Les n + 1 polynômes {L k )k=Q\ ...n son t linéairement indépendants. 

En effet soit oq, ai, ...cx n n +l réels tels que: 

a Q L (X) + aiLt(X) + . . . + a n L n (X) = 

La relation précédente devrait être vraie pour tout x G IR en particulier pour x = (^)i=o,i ! ... l n- 
pour X = Xi, i = 0, 1, .. . ,n ou a aiLi(xi) = donc Oui = , i = 0, 1,... ,n. La famille (£fc)fc=o,i,...,7i 
est une famille de n + 1 él'ements linéairement indépendantes, c'est donc une base de l'espace 
vectoriel des polynômes de degré n. 
Compte tenu de Qi et Q2, P s'écrit de manière unique comme combinaison linéaire des L k : 

P(X) = Y, 0kL k (X) 

comme P(xî) = fi i = 0,1, ... , n, on trouve donc que J3i = f%. 
Conclusion: 

P{X) = f L (X) + /iLi(X) + . . . + f n L n (X) 

Les L k sont appelés les polynômes de Lagrange élémentaires et P c'est le polynôme de Lagrange. 
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4.3.2 Exemple 

Nous allons construire le polynôme de Lagrange passant par les 4 points: 

A (2,i) ■ Ai (3, 2) , A 2 (-l,3) , A3 (4, 4) 
on pose: £0 = 2, #1 = 3, X2 = — 1 et xz = 4. /q = 1, f\ = 2, /2 = 3 et /3 = 4. 

Polynômes de Lagrange élémentaires: 

= (x-3)(.y-(-i))(x-4) = 1 _ _ _ 

w (2-3)(2-(-l))(2-4) 6 V A k " y ' 

= { X-2)(X-(-mX-4) = -1 _ _ _ 

w (3-2)(3-(-l))(3-4) 3 V A v m ' 

T I \ (X-2)(X-3)(X-4) -1 

L *) = (_i_ 2)( _i_3)(-l-4) = W (X " 2)(X " 3)(X " 4) 
i3(x) = (X-2)(X-Z)(X -(-!)) = ! _ _ _ 

iu (4-2)(4-3)(4- (-1)) W y A A v " 

Le polynôme de Lagrange est don donné par: 

P(X) = 1L (X) + 2L 1 {X)+3L 2 (X)+AL 3 (X) 

60 20 15 5 

4.3.3 Autre écriture du Polynôme de Lagrange 

1. Soit V(X) = (X - x )(X - xi) . . . (X - x„) = ngS(X - Xi), on montre facilement que: 

V'{x k ) = (x k - x )(x k - xi) ... (x k - Xk-i)(Xh - Xk+i)-.. (Xk -x n ) 

= (TV=2(Xk-Xi))^ k (4.3) 

Les polynômes de Lagrange élémentaires peuvent s'écrire comme: 

V(X) 



et par conséquent: 



Lk{X) ^V'(X)(X-x k ) 



p(x) - V fkV{x) 

P{X) ' Jinw-**) 



nX) k f vlf J k T _, (4.4) 



^ V'{X)(X - x k ) 

2. Variavles réduites 

Supposons que les x% sont équidistants: 

X{ = xq + ih , i = 0, 1,. . . , n , h = Xi+i — x% 
On introduit la variable réduite s de la façon suivante: 

X = Xq + sh 

Si X = Xq S = 

si X = x n s = n car x n — xq = nh 
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54 3. Interpolation et approximation polynômiale 

comme (X — Xi) = xq + sh — xq — ih = h(s — i) et (#& — x^) = xq + kh — xq — ih = h{k — i), L^ 
prend la forme suivante: 

L k (X) = L k (x + sh) = ngg {s ~ t} 

(fc - l)j^ k 

Le polynôme d'interpolation de Lagrange devient donc: 

k=0 l* W* 

4.3.4 Estimation de l'erreur 

Thôrème: Soit / une fonction de classe C n+1 [a,b], Si Xq } %±,. . . } x n £E [a,fr], -P le polynôme 
d'interpolation de Lagrange au points: (xq, f(%o)), (xi, f(xi)) : . . . , (x n: f{x n )) alors: 

Va- É [a, 6] 

/(*) " P(*) = ^f|^/ M («) (4-5) 

où a < miri(ar,Xo) < 9 < max(x, x n ) < b 

Preuve: Si X = Xi, f{x) — P(x) = donc l'équation 4.5 est vérifiée. 
Si x ^ Xi soit ip la fonction auxiliaire définie par: 

ij){t) = f{t) — P(i) — (t — Xo)(t — Xi) . . . (t — x n )(f){x) xest un paramètre 

<P( X ) = 7 û ï 7 \ 

(x -Xq){X- X1)...(X- X n ) 

on remarque que: iP(xq) = ij){x\) = . . . = ip{ x n) — i>{%) — •> donc ip est un polynôme ayant n + 2 
racines. 

Résultats préliminaires 

Si g [a, b] — > IR telle que g' existe: 

Si g possède (n + 2) zéros distincts alors g 1 possède au moins (n + 1) zéros distincts, 
généralisation du théorème de Rolle: Si g G C' n+1 '([a, b]) et si g possède au moins (n + 2) zéros 
alors g( n+1 ) a au moins un zéro. 

D'après le théorème de Rolle généralisé, la fonction ip( n+1 > admet au moins un zéro. Soit 9 ç]a, b[ 
tel que: ip i - n+1 '(9) = 0, On conclut donc que l'erreur est donnée par: 

= {x-x Q ){x- Xl ) (x-x n ) (n+1) 
w (n + 1)! w 

De plus si f( nJrl ) est continue sur [a, 6], elle atteint sa borne supérieure, on a la majoration de 
l'erreur suivante: 

TP ( \ \tt \ U( M ^ \{X-X { ))(x-X 1 )...(x-X n )\ r(n+l)f U 

e n{x) = f(x)-P(x) < maxf + '(as)) 

[n + 1)! xe[o,6] 

Exemple: Construire sur [—5, 5], les polynômes d'interpolation de Lagrange de la fonction / définie 
par f(x) = .._} 5 en prenant des abscisses d'interpolation équidistantes. 

On constate que l'on approche bien la courbe 1 2 - au centre de l'intervalle. Par contre, aux 
extrémités se produisent des oscillations parasites appelées effet de bord. Nous allons voir dans la 
section 6 que ses effets de bord disparaissent si l'on choisit comme points d'interpolation les points 
de Chebyshev. 
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Figure 4.2: Interpolation Poîynômiale de Lagrange. 



4.4 Interpolation poîynômiale de Newton 

4.4.1 Polynôme de Newton 

L'inconvénient des polynômes de Lagrange c'est qu'il n'y a pas de relation de récurrence entre eux. 
En d'autres termes, chaque polynône doit être construit individuellement et il n'y a pas moyen 
de calculer le polynôme de degré n à partir de celui de degré n — 1: les deux polynômes sont 
complètement indépendants. La méthode de Newton nous permettra de déduire P n à partir de 
Pn-1- Cette méthode est introduite récursivement de la façon suivante: 

Pi{x) =a + a 1 (x-x ) (4.6) 

Pi(x) = a + a\(x - x ) + a 2 (x - x a )(x - x\) (4.7) 

P 3 (:r) = a + ai(x - x ) + a 2 (x - x )(x - xi) + a 3 (x - x )(x - xi)(x - x 2 ) 



Pn(%) = a + ai(x - Xq) + a 2 (x - x )(x - xi) + a 3 (x - x )(x - xi)(x - x 2 ) 
+a 4 (x - x )(x - xi)(x - x 2 )(x — ses) + ... + 
o,n{x - xo) ■ • • {x - ïjv-i) (4.8) 

On remarque que Pn est obtenu à partir de Pn-\ en utilisant la relation: 

P N {x) = P N _i(x) + a N (x- x ){x- Xi)...(x- kjv_i) 

Le polynôme Pjy (4-8), de degré < N, est appelé polynôme de Newton au points xq 7 xi, . . .a;jv-i- 
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56 3. Interpolation et approximation polynômiale 

4.4.2 Exemple 

Etant donné les points xo = 1, X\ = 3, Xi = 4 et £3 = 4.5 et les coefficients ao — 5, a i — — 2, 
a 2 = 0.5, a 3 = -0.1 et a 4 = 0.003. Trouver Pi(l), P 2 (x), P 3 (x) et P 4 (x) et calculer P(2.5) pour 
i= 1,2,3,4. 
En utilisant les formules précédentes on a: 

P,(x) = 5-2(i-l) , P 2 (x) = 5-2(a;-l) + 0.5(a;-l)(x-3) 
P 3 (x) = 5 - 2(as - 1) + 0.5(2; - l)(x - 3) - 0.1(i - l)(x - 3)(x - 4) 
P 4 (x) = P 3 (s) + 0.003(i; - 1)(.t - 3)(x - 4) (x - 4.5) 

On trouve: Pi (2.5) = 2, 

P 2 (2.5) = Pi(2.5) + 0.5(1. 5)(-0.5) = 1.625, 

P 3 (2.5) = P 2 (2.5) - 0.1(1.5)(-0.5)(-1.5) = 1.5125 

et P 4 (2.5) = P 3 (2.5) + 0.003(1.5) (-0.5) (-1.5) (-2) = 1.50575 

4.4.3 Interpolation polynômiale de Newton 

Nous allons déterminer les coefficients a& pour tout les polynômes Pi, P 2 ,. . . , Pat qui approchent 
une fonction / donnée. La détermination des constantes a^ se fait à partir du fait que les polynômes 
de Newton Pj coï ncident avec la fonction / au points xq, xi, . . . , xy. f(xi) = Pj(xi) i = 0, 1,... ,j 

Pour le polynôme Pi, ao et ai sont déterminés à partir de 

P 1 (x ) = f(x ) et Pi{x l ) = î(x 1 ) 
On trouve: 

a = f(x„) (4.9) 

Xi - X Xi - x 

Les coefficients ao e t ai de P 2 sont les mêmes que ceux qui figurent dans Pi et qu'on vient de 
calculer (4.9,4.10). Le coeficient a 2 de P 2 est déterminé de la façon suivante: 

P2{X2) = f(x 2 ) = a + a 1 (x 2 - x ) + a 2 (x 2 - x )(x 2 - xi) 

a 2 est donc donné par: 

f(x 2 ) - a a - a 1 (x 2 - x ) 



11, 



(x 2 - x Q )(x 2 - Xi) 
1 J(x 2 ) - f( Xl ) /( a ;i)-/(x ) 1 



(4.11) 



(x 2 - x a ) X 2 - Xi Xi - Xo 

Cela nous amène a introduire le concept des différences divisées: 

Définition: Les différences divisées d'une fonction / sont définies comme suit: 

f[x k ] = f{x k ) 

t < 1 fl x k] ~ f[xk-l] 

f[Xk-i,x k \ = 

%k ~ Xk-l 

Cl 1 f[x k -\,Xk\ ~ /[xfc-2,Xfc_l] 

f[xk-2,Xk-i,Xk] = (4.12) 

%k - %k- 2 

,, , f[Xk-2,Xk-l,Xk\- f[Xk-3,%k-2,Xk-l] 

f[Xk-3,Xk- 2 ,x k -l,x k \ = 

Xk - Xfc_ 3 
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3.4 Poîynômiale de Chebyshev 

La formule de récurrence est donc définie par: 

f[x k _ i+1 ,...,x k ] - f[x k _j 



f[xk-i,Xk-i+l, ■ ■ ■ ,Xk] ■ 



>a;*-ij 



Xk ~ Xk- 
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(4.13) 



Théorème: Soit / 6 C n+1 ([a,b]). Si Xq, Xi, ..., x^ sont N + 1 points distincts de [a, b], alors il 
existe un unique polynôme P/v, dz£ polynôme de Newtom, de degré au plus égal à TV tel que: 

PN{Xi) = f(Xi) pour i = 0,l,...,iV 

P N (X) =a + ai(X - Xo) + ... a;v(X - :co)(^ - #1) .. . (X - Xn-i) 

où ai = f[xQ, x%, . . . ,Xi] pour i = 0,l,...,N 

de plus, pour tout x ÉE [a, b] il existe 9 ÉE [a, b] tel que: 

/(a:)=P !V (x) + B J v(x) 



w = (x ~* o)( y^y,- ( *~^V +1) ffl 



(4.14) 



Sjv est l'erreur de l'approximation de Newton. 

Exemple de construction du Polynôme de Newton 

Soit f(x) = x 4 — 2x 3 + X — X — 1, construisons les différences divisées basées sur les points xq = 1, 
#1 = 2, £2 — 3, X3 = 4, £4 = 5 et £5 = 6: 



Xi 


/k.l 


/[,] 


/[ , , ] 


/[ , , , ] 


/[.,,,] 


/[,.,.,] 


050 = 1 


-2 


- 


- 


- 




- 




- 


Xi = 2 


1 


3 


- 


- 




- 




- 


x 2 = 3 


32 


31 


14 


- 




- 




- 


x 3 = 4 


139 


107 


38 


8 




- 




- 


x 4 = 5 


394 


255 


74 


12 




1 




- 


x 5 = 6 


893 


499 


122 


1G 




1 








Le polynôme de Newton s'écrit donc comme suit: 

P 4 (x) = -2 + {x- 1)[3 + (x - 2) [14 + (x - 3) [8 + (x - 4)]]] 

4.5 Polynôme de Chebyshev 

4.5.1 Polynôme de Chebyshev 

La formule (4.14) montre que l'erreur de l'approximation poîynômiale de Newton est un produit 
de la dérivée (n + l) eme de /, évaluée à un point inconnu, avec l'expression (x — Xq) . . . (x — xn) 
qui ne dépend que de la subdivision (xq,Xi, . . . ,Xn). Il est alors intéressant de chercher, pour un 
n donné, la subdivision de [a, b] pour laquelle: 



soit minimal 
Le polynôme de Chebyshev apporte un élément de réponse à cette question. 



max \(x — Xq) . . . (x — xn)\ 



(4.15) 



Définition 

Pour n = 0, 1, 2, . . . et pour tout x G [—1, 1], les polynômes de Chebyshev sont définis par: 

T n (x) = cos(nArccosx). 
Il est facile de se convaincre que: Tq(x) = 1, T\{x) = x et T%{x) = 2x 2 — 1. 
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Propriétés des Polynômes de Chebyshev 

i) Les fonctions T n (x) satisfont la relation de récurrence suivante: 
Tn+i(x) = 2xT„(x) - r„_i(a;). 



(4.16) 



Par conséquent, T n (x) est un polynôme de degré n dont le coefficient de x n est 2" , c'est à 
dire, T n {x) = 2 n - 1 x n ... 

ii) \T n (x)\ < 1 pour x E [-1,1], T 2m (-x) = T 2m (x) et T 2m+ i(-x) = -T 2m+1 (x) 
iii) T„(cos(*f )) = (-1)*, pour k = 0, 1, . . . ,n 
iv) T„(cos( (2fc 2 + n 1)T )) = 0, pour fc = 0, 1, . . . , n - 1 

v) Les polynômes T n (x) sont orthogonaux par rapport à la fonction de poids -7= — y: 



/; 



(/.r 



VT- 



=T n (x)T m (x) 



n si n = m = 
| si n = 777, 7^ 
si n/m 



Preuve 

i) l'équation (4.16) est une conséquence immédiate de la relation: 

cos((?7 + l)#) + cos((tt,- 1)0) = 2cos(/9)cos(n6>). 

ii), iii) et iv) sont faciles à établir 

v) Pour montrer l'orthogonalité, on pose cos(#) = x et donc = Arccosjr: 



i: 



<i.r 



VI -x 2 ' 



%(x)T m (x) = / cos(n&) cos(mÔ)dO 
■h) 

1 f 7 " 

= - i [cos((n + m) 6) + cos((n — m 

2 ./o 

(n si n = m = 
| si n = m / 
Si 77 7^ 777. 

A l'aide de la relation (4.16) on montre facilement qu'on a le Tableau suivant: 

T (x 

Ti{x 

T 2 (x) = '_>./-' 

T 3 (x 

Ti(x 

T 6 (x) = 2 6 ~V 

T 6 (a;) = 2 fi - 1 x 6 

T 7 (x) = 2 7 ' 



■ 3x 

■ 8x 2 + 1 

■ 20a! 3 + 5x 

■ 48x 4 + 18x 2 - 1 

■ 112x 5 + 56x 3 - Ix 



(4.17) 



(4.18) 



Remarque: 

D'après la propriété i), T n est un polynôme dont le terme de plus haut degré x n a pour coefficient 
2 n . On appelle polynôme de Chebyshev normalisé à l'unité de T n (x) le polynôme T n (x) défini 
par T n (x) = 2 J- l 'Info) et on a: 

max T n (x) ■ 
a:e[-i,i] 



2 n-i 
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Théorème: 

Si E n désigne l'ensemble des polynômes de degré n et dont le coefficient de x 11 est 1, alors ponr 
tout polynôme p de E n on a: 

r = max T n (x) < max \p(x)\ (4.19) 

2"" 1 a, e [-i,i] "^£[-1,1] 

Preuve: 



Supposons qu'il existe un polynôme p de _E n tel que: max \p(x)\ < — — r et considérons le 



1 
xipï;i] '" v ~ " 2» 

polynôme rf(.x) défini par: d(#) = T n (x) — p(x). d est un polynôme de degré < n — 1 (le coefficient 
de x n est le même pour T n {x) et p{x)). 
La fonction polynômiale d{x) s'annulle au moins une fois dans chacun des intervalles fermés: 

,&7T (k + 1)71",-, 

Ixfc = cos( — ), Xk+ 1 = cos( )| , fc = 0, 1, . . . ,n — 1 

n n 

En effet: d'après iii) on a: d(xk) = ^=1 p( x k) et d(xk+i) = 2W _i p(#fc+i)î il est facile de 

se convaincre que d(#fc)d (#&.{. 1) < 0- Alors, <i(ar) possède n zéros dans [— 1, 1] et degd < rt— 1, par 
conséquent d(x) est identiquement nul donc: 

d(x) = , T n (x) = p(x) 

= max |T„(;r)| = max \p(x)\ < (4.20) 

ce qui est en contradiction avec l'hypothèse: max |p(aOI < r . 

1 JP a:e[-l,l]^ 2"" 1 

Le théorème précédent montre bien que: 

max |(.t — xq)(x — xi) . . . (x — x n )\ est minimal 
xe[-i,i] 

si et seulement si (x — Xq){x — x±) . . . (x — x n ) = T n (x), c.à.d, si 

,(2A; + l)7r N 
x k = cos( v , \ ) , fc = 0,L...,n 
v 2(n + l) ; ' ' 



les Xfc son t appelés les points de Chebyshev. Pour répondre à la question (4.15), il faut utiliser la 

a+b 
■2 



translation x — > ^— H — Fp£ qui envoie l'intervalle [— 1, 1] sur [a, 6]. On a donc le théorème suivant: 



Théorème: 

L'expression (4.15) est minimale parmi toutes les subdivisions (jpQ,aïi, . . . ,x n ) si et seulement si: 

a + b b- a ,(2fc + l)-7ï\ 
x k = + cos P , \ ) , k = 0,l,...,n. 

k 2 2 V 2(n + l) ' ' ' ' ' 

maXa-ei^Ka: - xq){x -xi)... [x - x n )\ = ^ 2n+1 — . 



(b-a) 



n+l 



\f(x)-Pc(x)\ < -^ ' , max |f ln+1 '(x)| (4.21) 

lJy ' ov /J — (n + l)!2 2 "+ 1 ïgI-1,1] 

Pc est le polynôme d'interpolation de Lagrange basée sur les points de Chebyshev. 

Exemple: Contruire sur [—1,1] les polynômes d'interpolation de Lagrange de la fonction f(x) = 
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Figure 4.3: Polynôme de Chebyshev. 



tttô .2 d'une part en prenant des abscisses d'interpolation équidistantes, d'autre part en choisissant 
les points de Chebyshev. 

Cet exemple est dû a Runge (appelé phénomène de Runge) met bien en évidence les effets 
de bord. Par contre, si on utilise les points de Chebyshev, on constate que les effets de bord 
disparaissent. 
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4.6 Exercices 

Exercice 1. Si f(x) = e x , montrer que n = 9 est le plus petit entier telle que l'erreur \Eg(x)\ = 
\f{x) - P 9 {x}\ < 0.0000005 pour x G [-1. 1] 

Exercice 2. 

a. En utilisant la série géométrique j^ = 1 — t 2 + t A — t® -\- . . . \t\ < 1 

après intégration montrer que : ArcTan(.T) = x — ^r + ^ — y- + • • * 

b. Montrer que: n = 2^/3(1 -Ç + Ç-^r + Ç----) 

c. Utiliser b. pour calculer la valeur approchée de n à 10 _s . 

Exercice 3. Soit f(x) = ln(l + x). 

a. Montrer que /(*) = (-l)^" 1 )^ - 1)!/(1 + x) k 

b. Ecrire le développement de Taylor de / à l'ordre n. 

c. Calculer P 3 (0.5), P 6 (0.5) et P 9 (0.5). Comparer à /(0.5) 

d. Montrer que pour x ÉE [0,0.5], la majoration de l'erreur est donnée par: E$(%) < 0.00009765 

Exercice 4. 

1. Ecrire le polynôme d'interpolation de Lagrange P de la fonction f(x) = ln^ pour les 4 abscisses 
d'interpolation suivantes: 1, 2, 3, 4. 

2. Calculer P(tt). Comparer cette valeur à la valeur exacte In j. Donner une majoration de l'erreur. 

3. Calculer P'(4) et comparer cette valeur à /'(4). 

Exercice 5. Soit g G C +1 [a, 6], supposons qu'il existe deux points x 7 xq ÉE [a, b] tels que: 
g(x) = g(z Q ) = g'(xo) = g"(x Q ) = ...= g^ N) {x ) = ° 

a. Montrer qu'il existe un point c entre x et xq tel que: g( N+1 >(c) = 

b. Soient / 6 C'' iV+1 '[a,&], Pjv le développement de Taylor au voisinage de xq à l'ordre TV de / et 
Rn le reste du développement tels que: 

f(x) = Pjv(#) + Rn(x,xq,N : Ô), est compris entre x et xq. 



Ecrire le polynôme P/v- 
c. Démontrer que 



R N (x,x Q ,N,8) = f( N + 1 \ê) { - X Xo)mi 



{N + iy. 

Exercice 6. 

L'objet de l'exercice est le calcul approché de l'intégrale d'une fonction continue / de [a, b] dans IR 
1. / désigne une fonction de classe C 2 de [a, b] dans IR. On pose pour k = 0, 1, . . . ,n: %k = a+&^p. 
Soit P le polynôme de degré au plus égal à 1 tel que P(a) = /(a), P(b) = f(b). 

a. Donner la valeur de l'intégrale de P sur [a, b]. On note / cette valeur. 

b. Par application du théorème de Rolle à une fonction auxiliaire convenablement choisie, prouver 
que pour tout x G [a, b], il existe c £]a, b[ tel que: 

f( x )-P( x ) = -t±^( x -a){x-b) 

et en déduire que si Ma est un majorant de |/"| sur [a, b] alors: 

r b Ah 



/0 TUT 

f( x )dx-I\<^{b-af 
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2. Par application de cette méthode sur chaque segment [$&— i, Xk] (1 £ k < ^) établir la majoration 
suivante: 

1/ /(x)<fa-/„|<-^(6- a ) 3 

b — a . f(à) ~^K .. , /(&)-. 

avec I n = (^+ X] /fe) + ^) 

ri z fc=1 ^ 

Exercice 7. Soient (ûi)i=o.i,...,n une subdivision régulière de [a, 6] de pas h et / € C n+ [a, &]. 
Considérons les points j4ï(Oz, /(ai))z=o,i....,n de IR , on désire construire un polynôme P n (X) de 
degré au plus égal à n passant par les points Ai{di, /(ai))i=o,l,...,n- 

a. Démontrer que les polynômes (Ij(X))j = q^ n de degré n tels que: L'(ttj) = <5jj i = 0, 1, . . . , n 

sont de la forme: 

b. En utilisant le changement X — qq = th montrer que: 



n?=o(*-0 , 
n[L 0' - 



W = ;t° _ 1^- ^O' 1 ' 



Si .t G [a, b] quel est le domaine de variation de t? 

c. Montrer que la famille (Jq(X),Ii(X), . . . ,l n (X)) est une base de l'ensemble des polynômes de 
degrés < n. 

d. Démontrer que le polynôme de Lagrange P n est de la forme: 

P n (X) = XÏ =0 f( ai )k(X) 

e. En appliquant le théorème de Rolle à une fonction auxiliaire convenablement choisit montrer 
que: 

P n (x)-f(x) = {x ~ Qn)(X (t ; + ° 1 1 ) )! - - (X ~ a "V " +1) ffl = E n (X) , ee]a,b[ 



et en déduire que 



{n + lj! 



f. Montrer que si M n+ i est un majorant de f( n+1 > sur [a, b] alors 

(t - 2) ■ ■ 
(ra + 1)! 



P n{X) - f (x)\ = \E n{X )\ < r» t(t - 1)( *- 2] ';;- (f -" ) M„ +1 



Application: Montrer que: 

, / x. h 2 M 2 

\Ei{x)\ < —z— xe [a 0) «i] 

o 

, s h 3 AU 
\E 2 {x)\<—j- iE|o,,flJ 

, ,. fr 4 i\/ 4 
l-^aWl < — ^p x e [a ,a 3 J 

Exercice 8. Soient lk{x)(k=0,l,2) l es polynômes de Lagrange élémentaires au points xq, x\ et x%. 
On définit g par g(i) = Zo(^) + h( x ) + fe(ï) — 1 
a. Montrer que g est un polynôme de degré < 2. 
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b. Montrer que g(Xk) = pour k = 0.1.2 

c. Montrer que g(x) = pour tout x. 

d. généraliser le résultat au cas de n racines Xq, X\. ■ ■ X%. 

Exercice 9. Soient lk( x )(k=o,l,...,n) l es polynômes de Lagrange élémentaires au points Xq, Xi, ..., 
x n . On définit g par g(x) = Yli=Q x iLi{x), montrer que pour tout x G IR, g(x) = x. 

Exercice 10 : a. Soit / G (7 3 [a, 6], xq,xi,X2 trois points de [a, b]. 

Soit Pi(X) = an + a,\{X — Xq) le polynôme d'interpolation de Newton du premier degré passant 
par Aq(xq, /(xq)) et A\(x\^ f(x\)). Calculer clq et ai en fonction de xq, X\ et f(x$) : /(xi). 
b. Soit P2{X) = P\(X) + a^(X — xq)(X — x\) le polynôme d'interpolation de Newton du second 
ordre passant par Aq(xq, /(xq)), Ai(xi, f(x±)) et A^ixi-, f(x%)). Montrer que: 

,, , /(xi) - f(x ) 



a-2 ■ 



Xl - X 
1 f( X2 )-f( Xl ) f( Xl )-f( XQ ) 



(x 2 - Xq) x 2 -Xl Xl- Xq 



c. En utilisant une fonction auxiliaire convenablement choisie, montrer que pour tout x G [a, b] on 
a: 

f(x) = P 2 (x) + (x - x )(x - x x )(x - ^Y—^- 1 - = ft(x) + E 2 (x) c e]a, b[ 

d. Supposant que xq, x\ et X2 sont equidistants, Exprimer B^Çx) en terme de la variable réduite s 

(x = Xq + sk) et h = Xi — Xq. 

e. Soit M3 un majorant de f^' sur ]a, o[, montrer que 

1^2(^)1 < 



9^3 

Exercice 11 : On rappelle la méthode de Lagrange. qui à partir des 2 valeurs initiales Xq et Xi, 
permet de générer une itération afin de résoudre l'équation f(x) = 0, soit: 

x n+l — *ïî J \ x n 



f( x n)~ f{ x n-l) 

a. Interpréter géométriquement cette méthode en montrant que x n+ j est l'abcisse du point 
d'intersection de la droite sécante passant par les points (2^-1, f(x n -i)) G Cf et (x n , f(x n )) G Ct 
avec l'axe des x. Faire une figure illustrant ce résultat dans un cas partculier. (Cf est la graphe de 

/) 

b. On peut également interpréter ce résultat en terme d'interpolation inverse. Ecrire le polynôme 
d'interpolation de Lagrange Q basé sur les points (f(x n —i),x n -i) G Cr-i et 

(f(x n ),x n ) G Cf-i- Calculer Q(0). Que retrouve-t-on? (C*-i est le graphe de / -1 ). 

c. Proposer une amélioration de la méthode précédente utilisant 3 valeurs successives de l'itération 
et donc un polynôme de degré 2. Calculer la valeur obtenue alors pour x n +i en fonction de x n , 
X n -i , x n _ 2 et f(x n ), f(x n -i), f(x n _ 2 )- 

Exercice 12. (Interpolation d'Hermite) Soit / une fonction numérique définie sur [a, bb] et /' 
sa dérivée. (iCj)j— o,,„,n une subdivision de [a, b]. Soient (fi = /(#$)) une suite de valeurs de / aux 

points (iCi) î - = o ! ... ! n et (yi = f'(xi)) une suite de valeurs de /' aux points {Xj)i=Q n . On rappelle que 

le polynôme de Lagrange d'indice k est donné par: 

T , v . (X- x„)(X - Xl )...(X- x k -i)[X - x k+1 ) ...(X- x n ) 



{x k - x )(x k -Xl)... (x k - x k _i)(x k - x k+1 ) ...{x k - x n ) 
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1.) Calculer le polynôme d'interpolation de Lagrange de f(x) = e~ x aux points 0, 1, 2. 
2.) On suppose qu'il existe un polynôme P de degré au plus 2n + 1 tel que: 

P(x l ) = fi et P'{xi) = fi pour i = 0, . . . , n (4.22) 

Montrer l'unicité d'un tel polynôme. 

3.) Pour fc = 0, . . . ,71, on considère les polynômes: 

iP k (X) = (X - x k )L 2 k {X) 

MX) = [1 - 2(X - x k ) dI ^l(x k )}Ll{x) 

et 

n 

Pm = J2lhMX) + ytMX)} 

3a.) Calculer iftk(%i)i ^fcO*^)' (fikixi), ^(^i) pour < i < n et en déduire que pour < i < n 

Pfe) = h et P'fc) = ./;' 

3b.) Quel est le degré de ip^, 0^ et P. Le polynôme ainsi construit est appelé le polynôme 

d'interpolation d'Hermite de / aux points (xi). 

4.) On supose que / € C 2n+2 [a, 6], montrer qu'il existe Ç € [a, b] tel que: 

/m - p(x) = ^-^m^.. { x- Xn r fi2n+2){0 

5.) Application: Calculer le polynôme d'interpolation d'Hermite de f(x) = e~ x aux points 0, 1, 
2. 
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